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Abstract 

The simplest statistical-mechanical model of crystalline formation (or alloy 
formation) that includes electronic degrees of freedom is solved exactly in 
the limit of large spatial dimensions and infinite interaction strength. The 
solutions contain both second-order phase transitions and first-order phase 
transitions (that involve phase-separation or segregation) which are likely to 
illustrate the basic physics behind the static charge-stripe ordering in cuprate 
systems. In addition, we find the spinodal-decomposition temperature satis- 
fies an approximate scaling law. 
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I. INTRODUCTION 



The most fundamental problem in solid-state physics is to understand why elements (and 
most compounds) crystallize in ordered periodic structures, for this forms the basis of all of 
solid-state physics. While it is well known that the driving principle behind this ordering is a 
lowering of the ground-state energy of the material, and there has been significant progress 
with ab initio methods to predict the ground-state properties of these ordered phases in 
real materials, there still are no exactly solvable models for crystal formation that describe 
the statistical-mechanical mechanism behind the ordering of the electrons and ions on a 
periodic lattice. Furthermore, it is not understood what the physical mechanisms are that 
are necessary for creating a crystallized state. This crystallization problem is ubiquitous; it 
also describes the statistical mechanics behind binary-alloy formation or phase separation 
since the two problems can be mapped onto each other (as described below), and it may 
also describe the physics behind charge-stripe formation in the cuprates. 

It may sound surprising that no solvable statistical-mechanical model for crystallization 
exists, since a statistical-mechanical model for magnetic order has been known ever since 
Onsager solved the two-dimensional Ising mo Onsager's solution produced a paradigm 
for understanding phase transitions in many different physical systems and provided a text- 
book example of much of the theory behind modern critical phenomena. In fact, Lee and 
Yang§ modified the Ising model to consider the magnetic order in an external magnetic field, 
and mapped the problem onto a lattice gas, where the up spins denoted sites occupied by 
ions, and the down spins denoted empty sites. Onsager's method of solution does not extend 
to the case of a finite magnetic field, so no exact results are known for the lattice gas, except 
in the case where the number of ions equals one half the number of lattice sites, which corre- 
sponds to the zero-field case. These models of crystallization neglect the electronic degrees 
of freedom of the valence electrons, and hence are not directly applicable to real materials 
such as metals and alloys. 

It turns out that the Ising model, and many other models for magnetism, simplify when 
they are examined in high dimensions. In fact, the Ising model is solved by a static mean 
field theory in four and higher dimensions. A similar situation is expected for electronic 
problems, except they remain nontrivial even in the infinite-dimensional limitiS. Metzner 
and Vollhardt showed that the electronic problem requires a dynamical mean-field theory 
for its solution in infinite dimensions. Furthermore, a wide range of evidence indicates that 
this dynamical mean-field theory provides a quantitative approximation to the solutions of 
correlated electron problems in three dimensions (at least if one is not too close to a critical 
point). In fact, it is precisely the nonuniversal properties (such as a transition temperature) 
that the dynamical mean-field theory determines accurately, and it's solution provides a 
wealth of information on the qualitative behavior of the model studied. We employ the 
dynamical mean-field theory here to produce an exact solution of the crystallization problem 
which includes the electronic degrees of freedom. 

The simplest model that can describe crystallization and include electronic degrees of 
freedom is the spinless Falicov-Kimball mo del which consists of two kinds of particles: lo- 
calized ions and itinerant (spinless) electrons. The localized ions (wi = or 1) occupy sites 
on a lattice in real space with an energy E, and the electrons can hop (with a hopping inte- 
gral — t*/[2\/d\) between neighboring lattice sites. In addition, there is a screened Coulomb 
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interaction U between electrons and ions that occupy the same lattice site. Since the elec- 
trons do not interact with each other, the "spin" degree of freedom is unimportant, and is 
neglected. The Hamiltonian is 

t* 

H = ~—y= ^ C Uj + E Y, W i + U Y, C i C i W i> (!) 
zy a <jj> i i 

with cj (q) the creation (annihilation) operator for electrons at site i, and Wi denoting the 
ion occupancy at site i. We use t* = 1 as the energy scale. 

The Falicov-Kimball model can be viewed as a simplified approximation to a real material 
in a variety of ways. If the material has a single valence electron, and only one electronic 
band lies near the Fermi level, then the crystallization problem would correspond to the 
case where the electron and ion concentrations (p e and pi) are the same (which is called 
the neutral case), since one electron is donated by each ion. If, instead, there are many 
bands near the Fermi level, then one can map the combined bands into a single "effective" 
band which will have an electron filling determined by the average filling of the electrons in 
the most important band. In this case, each ion may donate only a fraction of an electron 
to the crystal, because the rest of the electron goes into other hybridized bands that lie 
close to the Fermi level. Hence, one may find it useful to also consider nonneutral cases for 
the crystallization problem, where the electron and ion concentrations are not equal. This 
model can also be mapped onto the binary alloy problem, where a site occupied by an ion 
is mapped to a site occupied by an A ion and a site unoccupied by an ion is mapped to a 
site occupied by a B ion, and the screened Coulomb interaction is mapped to the difference 
in the site energies for electrons on an A ion versus ona5 ion. 

As it stands, the Falicov-Kimball model doesn't appear to be a many-body problem at all, 
since the ions are localized and do not move, which implies that the quantum-mechanical 
problem for the electrons can be solved by diagonalizing a single-particle problem of an 
electron moving in the potential determined by the given configuration of the ions {wi}. 
The many-body problem aspects enter by taking an annealed average over all possible ion 
configurations with the chosen ion concentration. This produces long-range interactions 
between the ions, that can cause them to order or phase separate at low temperatures. 

Much is already known about the physics of the Falicov-Kimball model (as reviewed 
by Gruber and Macrisil). In the neutral case where each particle concentration equals 1/2, 
Lieb and Kennedy! and Brandt and Schmidt! proved that the system always orders in an 
alternating "chessboard" phase at a finite transition temperature in all dimensions greater 
than 1. This ordered phase can be interpreted as the transition from a high-temperature 
homogeneous (liquid/gas) phase to a low-temperature ordered (solid) phase. The appearance 
of a low-temperature ordered phase follows as a consequence of the Pauli exclusion principle, 
since Lieb and Kennedy also showed that if the itinerant particles were Bosons instead of 
electrons, they would clump together and not form a periodically ordered ground state. 

The Falicov-Kimball model is expected to be in the same universality class as the Ising 
model, but, because of the electronic degrees of freedom, one needs to solve the full statistical 
model to determine the "effective magnetic exchange parameters" between different lattice 
sites. The parameters can be extracted in a systematic expansion if the electronic kinetic 
energy (the hopping term) is taken as a perturbationJHliiJ but such an analysis is only 
valid in the strong-coupling regime, and rapidly becomes problematic. It is precisely this 
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complication that has frustrated attempts at finding an exact solution to the crystallization 
problem when electronic degrees of freedom are introduced. 

The one-dimensional limit of the Falicov-Kimball model has also been extensively stud- 
ied. Here there are no finite-temperature phase transitions, but the system can have phase 
transitions in the ground state. The first attempt at studying the one-dimensional Falicov- 
Kimball model proceeded along the lines of ab initio band-structure calculations for real 
materials — a small number of candidate ion configurations were chosen for the ground state, 
and a restricted phase diagram was determined for all structures within the subsetEl. The 
numerical solutions produced two conjectures: the first was a result for the case where 
p e 7^ 1 — pi, which stated that if the screened Coulomb interaction U was large enough, then 
the system would segregate into an empty lattice (with no ions and all the electrons), and 
a full lattice (with all the ions and no electrons). The second was a generalization of the 
Peierls instability, which says that in the small U limit the system will order in such a fashion 
that the ions produce a band structure that has a maximal gap at the Fermi level. This 
first conjecture (the segregation principle) was later proven to be true by Lembergerll! while 
the second conjecture was shown to be false if the electron concentration was sufficiently far 
from half-filling. In that case, the system would phase separate between the empty lattice, 
and an optimally chosen ion structure that had the Fermi level lying in the gapEj. 

The other limit that has been extensively studied is the large- dimensional limit where 
Brandt and MielschEl provided the solution of the transition temperature as a function 
of U for the half-filled symmetric case. Their solution involves solving a coupled set of 
transcendental equations which display first and second-order phase transitions. Freericks0 
later showed that the model (on a hypercubic lattice) also displayed incommensurate order 
and segregation. 

There are two kinds of lattices that are usually investigated in the large coordination- 
number limit: the hypercubic lattice, which is the generalization of the cubic lattice to large 
dimensions; and the Bethe lattice, which is a thermodynamic limit of the Cayley tree when 
the number of nearest neighbors becomes large. The noninteracting band structure for the 
hypercubic lattice produces a density of states that is a Gaussian [p#(e) = exp(—e 2 )/y/7r\, 
while on the Bethe lattice the density of states is Wigner's semicircle [ps(e) = a/4 — e 2 / (2tt)]. 
The hypercubic density of states has an infinite bandwidth, but most of the weight lies within 
a range of ±2 about the origin. The Bethe lattice density of states has the same behavior 
as a three-dimensional system at the band edge (square-root behavior) but has no van Hove 
singularities in the interior of the band. Because both density of states are nontrivial, the 
many-body problem maintains much of it's rich behavior that arises from the competition 
between kinetic-energy effects and interaction-energy effects. In particular, the Falicov- 
Kimball model continues to have phase transitions in the large coordination number limit, 
but the transitions have mean-field theory exponents. 

In this contribution, we examine what happens in the case when the Coulomb interaction 
becomes infinite U — > oo (the attractive case is equivalent to this case through a particle-hole 
transformation of the electrons, which carries p e — > 1 — p e ). In this case, the electrons avoid 
the sites of the lattice occupied by the ions, so the electron concentration varies from zero 
up to 1 — pi. We investigated the non-unit-density cases, where the electron concentration 
was restricted to < p e < 1 — p,. In Section II the formalism and results for calculations on 
the Bethe lattice are presented. In Section III, results for the hypercubic lattice are given 
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and in Section IV we present our conclusions. 



II. FORMALISM AND RESULTS FOR THE BETHE LATTICE 

In the thermodynamic limit, the local lattice Green's function is defined to be 

G n = G{iu n ) = -J q dre ^ < ^h-^n) > > ( 2 ) 

where iu n = iirT(2n + 1) is the Fermionic Matsubara frequency, (3 = 1/T is the inverse 
temperature, p is the electron chemical potential, and T T denotes r-ordering. The angle 
brackets in Eq. (||) denote the sum over ionic configurations. The local Green's function is 
determined by mapping onto an atomic problem in a time-dependent field, with the following 
action 

S at = [ dr f dT'c\r)G G l (T - t')c(t') + U [ drJMcMw + Ew, (3) 
Jo Jo Jo 

where w — 0, 1 is the ion number for the atomic site and Gq 1 is the mean-field or effective- 
medium Green's function, which is determined self- consistently (as described below). The 
atomic Green's function, with the action in Eq. @, is computed to be 

with pi the average ion density < w >. On the other hand, the local lattice Green's function 
satisfies 



oc 



0„ = I de . P{€) , (5) 

iu n + n - h n - e 



oo 



where p(e) is the noninteracting density of states for the infinite lattice and E n is the self- 
energy. The self-consistency relation is that the self-energy S n in Eq. ([5]) must coincide with 
the self-energy of the atomic problem, i. e. 

Y,{iu n ) = G \iuj n )-G- 1 . (6) 

Equations and @ constitute the mean- field theory for homogeneous phases. In the 

limit d — > oo Eq. is an exact equation for the lattice problem. We note that for periodic 
phases, if they exist, one needs to replace the atomic problem by a more complicated many- 
site problemEJ. 

These equations are complicated to solve analytically but a simplification occurs for 
< p e < 1 — Pi in the limit U — > oo. Indeed, when U is large the spectrum of the 
Hamiltonian consists of two bands separated by a gap of order U for electron fillings that 
satisfy < p e < 1 — pi. In this case the chemical potential lies within the lower band, 
so that p is 0(1). We note that G is a function of p, and therefore for any finite p, 
l/[Go(iuj n )U] — ► as U — > oo. Then Eq. (|J) becomes 

G n = (1 - pi)G (iuj n ), (7) 
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and substituting Eq. (^) into Eq. @ and solving for the self energy, then yields 



= (8) 

for the relation between the local self energy and the Green's function. Hence, in the limits 
U — > oo and d — > oo the equations for the homogeneous phase reduce to Eqs. (|) and (|8|). 

In the case of the Bethe lattice, Pb(c) = V4 — e 2 / (27r), for — 2 < e < 2, so that the 
integral in Eq. (|5|) can be performed analytically 

icj n + //-£„ 1 /- ■ -— - 
G„ = - v /(^ n + /i-S n ) i -4. (9) 

Substituting the result from Eq. (§) into Eq. (|9D and solving for G n yields the exact result 
for the interacting Green's function in the strongly correlated limit 

iuj n + jj, 



G n = - -y/{iu n + /i) 2 - 4(1 - Pi ), (10) 

where the phase of the square root is chosen so that the Green's function has the correct sign 
to it's imaginary part. This form is identical to that of a noninteracting Green's function 
[Eq. (H) with E n = 0], with a bandwidth narrowed from 4 to 4a/1 — pi, and containing a 
spectral weight of 1 — pi (since the remaining spectral weight is shifted to infinite energies). 
This is easiest seen from the interacting density of states, which satisfies^ 

Pb%) = f^ftK. (ii) 

Note that in the infinite- interaction-strength limit, we have an analytic form for the Green's 
functions, and do not need to iteratively solve transcendental equations as is normally done 
in the finite-?/ caseS Furthermore, even though the Green's function has the same form as 
a noninteracting Green's function, the self-energy is nontrivial and does not correspond to 
a Fermi liquid! 

This form for the Green's function fits a rather simple physical picture. The electron 
avoids sites occupied by an ion when U — » oo, so the number of available sites is reduced by 
the fraction 1 — pi. This means, on average, the number of nearest neighbors is reduced by 
the same factor, which reduces the bandwidth by ^1 — p^. The total spectral weight is also 
reduced from 1 to l — pi, because the upper band (with p, t states) is located at infinite energy. 
What is surprising is that this "hand-waving" argument is exact for the Bethe lattice (we 
will see below it is a good approximation for the hypercubic lattice, but is not exact). 

The interacting density of states is temperature-independentl 2 ^ in the local approxima- 
tion, which means that we can examine the ground state at T = to see if the system phase 
separates, or if the homogeneous phase is lowest in energy. The ground state energy for an 
ion concentration p t and an electron concentration p e is 



E(p e , Pl )= I dep^(e)e, (12) 
with p the chemical potential defined by 



6 



dept{e) 



(13) 



and p % B l the interacting density of states. Substituting in the exact result from Eq. ( |TTD 
yields 



E(pe,Pi 



-f (1-P*) 3/2 

Sir 



1 3/2 



4(1 -Pi) 



(14) 



and 



Pe 



1 - Pi 



71 



COS 



-1 



-p 



p 



2^/T^p-) 2^T=7 l \ 4(1 -p 



P 1 



(15) 



Using Eqs. (|T4| ) and flT5|), we will show that the mixture of the state with no ions and an 
electron filling p e j (1 — p^ with the state with all ions and no electrons has a lower energy 
than the homogeneous state, i. e. 



E(p e ,p t ) > {l-fH)E( 



Pe 



1 - Pi 



,0) + Pi E{0,l). 



(16) 



Moreover, from Eq. ([U]) we will deduce that the mixture corresponding to the right hand 
side of Eq. (|T^) has lower energy than any other mixture between homogeneous states. In 
other words, 



aE{p' e ,p[) + (1 - a)E{plp'l) > (1 - Pl )E{ 



Pe 



1 - Pi 



0)+ Pi E(0,l), 



(17) 



where < a < 1 and p e = ap' e + (1 — a)p", pi = ap\ + (1 — a)p", < p' e < 1 — p[, and 
< p" < 1 — p". To obtain Eq. fll6|), we first notice that E(0, 1) = and that the chemical 
potential p corresponding to an electron filling of p e /(l — Pi) and an ion filling of zero is 
p = pj (2a/1 — pi), as can be seen from Eq. (|15|). Therefore, Eq. (|14D yields 



1 - Pi )E(-^-, 0) + Pi E(0, 1) = _ p . 

1 — Pi 07T 



p 



4(1 - Pi 



3/2 



:E(p e ,Pi) < E(p e ,Pi), 

(18) 



which proves Eq. (|IE|) . The proof of Eq. ( |TTD relies on an application of Eq. 

®E(p'e,Pi) + (1 - &)E{p'e,p") > 



a 



(i-p:)^( t A t,o) + p ^(o,i) 

1 - Pi 



p 



+ (l-a) 



l-pim-^-^ + p'lE^l) 

* Pi 



The right hand side of Eq. ( ]T5| ) is equal to 

;i-a)(l-tf) 



(1 - Pi) 



^1^(^,0) 
1 - Pi 1 - Pi 



1 - Pi 



p" 



(19) 



1-Pi 



+ Pi E(0,l). 



(20) 



On the other hand, E(p e , 0) is a convex function of p e , so the term inside the brackets 
in Eq. ( |20"1) is greater than E(p e /[1 — pi],0), which yields Eq. (|i~6"[) . We remark that the 
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convexity of E(p e ,0) is obvious from the fact that the free electron system cannot phase 
separate. Formally, it can be seen as follows: differentiating Eqs. ([14]) and fll5|) with respect 
to p e gives E'(p e ,0) = ppB(p)dp/dp e and 1 = pB(p)dp/dp e . Thus E'(p e ,0) = p and 
E"( Pe ,0)=dp/dp e = l/p B (p)>0. 

Our interest now is to determine the finite-temperature phase diagram of the infinite-t/ 
Falicov-Kimball model since we know the system always phase separates at low temperature 
(although we have not yet ruled out the possibility of charge-density-wave phases being 
lower in energy than the phase-separated ground state). The first step is to evaluate the 
conduction electron charge-density-wave susceptibility. It is often stated that the Bethe 
lattice can only support antiferromagnetic or uniform order — no incommensurate or other 
"periodic" phases can exist. But this statement has never been proven, and recent work 
has shown it to be false0 by a counterexample of a period-three phase stabilized on the 
infinite-dimensional Bethe lattice at zero temperature. The momentum dependence enters 
the dressed susceptibility only through the momentum dependence of the bare susceptibility 
because the vertex function is local in the infinite-dimensional limit. This allows us to simply 
take the U — > oo of the Brandt-Mielsch result^, which gives 



oo 



n=~ oo 



with X being the parameter that determines the modulation of the charge-density-wave over 
the Bethe lattice and with Xn(X) the corresponding bare susceptibility. We do not provide 
the general formula for all possible charge- density waves here. Rather, we present the three 
simplifying cases for the susceptibility on the Bethe lattice: (i) the local susceptibility, where 
Xn(local) = —G 2 n \ (ii) the [X = —1) "antiferromagnetic" susceptibility, where 

x° n (-i) = - . f n y ; (22) 

iuj n + p - S n 
and (iii) the (X = 1) uniform susceptibility, where 

XniX) = = 1, 23 



The local susceptibility never has a transition, because the numerator of Eq. (21) vanishes. 
The condition for an "antiferromagnetic" charge density wave becomes 

OO /~i 

l = ftE^ (24) 

after substituting in the infinite-?/ form for the self-energy, and using the quadratic equation 
G\ — (iu n +p)G n +l— pi = that the interacting Green's function satisfies. Now, substituting 
the integral form for G n 

G n = (l- p.) r de , PB(e) ^= , (25) 

J-oo lUJ n + p - Vl - Pi€ 

into Eq. ( p4|) and performing the sum over Matsubara frequencies yields the final integral 
form for T r 



S 



2T 7-2^7 2 cosh 2 -tanhf tanhf )' 1 ' 

(see Appendix A). But this integrand is positive for all z, so the right hand side is always 
less than zero, and there is no "antiferromagnetic" T c . The staggered charge-density-wave 
order has been found near half- filling Pi = p e = 1/2 when the lowest-order exchange for 
finite-f/ is included,^ but can only occur at T = and pi = p e = 1/2 when U = oo. 

The uniform susceptibility case is analyzed as follows: First the uniform susceptibility 
from Eq. (^) is substituted into Eq. (|2l|), and the square root is eliminated by using the 
exact form for the interacting Green's function in Eq. (|9|). Next, the self energy is replaced 
by its exact form from Eq. (|||), and the quadratic equation for G n is used to simplify the T c 
equation to 



Pi Yl 



1 — Pi 

1 + 



(iu n + p)G n - 2(1 - pi 



(27) 



Now the interacting form for G n from Eq. fllPl) is substituted into Eq. fl2T|) and the results 
simplified to yield 

_ - ( luJn + fl)Gn -2(l- Pi ) 

"Jko (-n + /i) 2 - 4(1 - Pi ) ■ {2S) 

The final step is to substitute in the integral form for G n from Eq. (^) and perform the 
summation over Matsubara frequencies (see Appendix A). After making a trigonometric 
substitution, the transcendental equation for T c becomes 

1 = Pi ^ n ~ Pi ^ 7r ^cosgtanh^(2 v /l - Pi cos 6 - fj). (29) 

We do not discuss any of the other periodic cases here, because the numerics involved is 
cumbersome. But, we expect the Bethe lattice to have similar behavior as the hypercubic 
lattice, where the transition always went into the uniform charge-density wave, signifying a 
phase separation transition. Details of the other periodic phases will be reported in a future 
publication. 

The results for the transition temperature for the uniform charge- density wave are pre- 
sented in Figure [TJ(a). We choose nine different ion concentrations ranging from 0.1 to 0.9 
in steps of 0.1. The electron density then varies from to 1 — p^ for each case. As can 
be seen in the figure, the maximal transition temperature is about 0.12t* and it occurs at 
half-filling of the lower band p e = (1 — Pi)/2 with pi m 0.65 (coincident ally, this maximal 
transition temperature is nearly identical to the maximal T c to charge- density- wave order 
at p e = pi = 1/2 when evaluated as a function of the interaction strength U). Since T c 1, 
we expand Eq. (p9"D for small T by replacing the tanh x by sgnx to find 



J-c ~ 



TT 



N 1 " W^TY (30) 



Since the chemical potential will scale with y/1 — pi for the same relative electron filling 



in the lower band [p e /(l — Pi)], as shown in eq. (|T5|) , this form motivates a scaling plot of 
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T c j (piy/1 — Pi) versus p e /(l — Pi), which appears in Figure 0(b). As can be seen there, the 
data nearly collapse on top of each other for T c which is usually a nonuniversal quantity. In 
fact, the variation in T c is less than 10% for all different cases. 

The susceptibility analysis shows that the system orders in a uniform charge-density- 
wave, which indicates that the system will phase separate (or segregate) into two regions, 
one with a higher concentration of electrons and one with a lower concentration (as we 
already showed at T = 0). Such a phase separation is usually associated with a first-order 
phase transition, rather than a second-order transition. Hence, it is important to perform a 
Maxwell construction of the free energy that includes mixtures of two states with different 
electron and ion concentrations such that p e = ap' e + (1 — a)p", p% = ap^ + (1 — a)p", and 
that the free energy of the mixture F(p' e ,p£; p",p") = aF(p' e ,p' i ) + (1 — a)F(p", p") is lower 
in energy than the pure-phase free energy F(p e ,pi). The second-order phase transition 
is the spinodal-decomposition temperature, below which the free energy becomes locally 
unstable in the region of (p e , pi); in most cases the global free energy is minimized by the 
Maxwell construction at a temperature above this spinodal-decomposition temperature. The 
spinodal-decomposition temperature marks the lowest temperature that the system can be 
supercooled to before it must undergo a phase transition. 

We can calculate the free energy F(p e ,pi) for a homogeneous phase with electron fill- 
ing p e and ion concentration pi in two equivalent ways. The first method is from Brandt 
and MielschB which expresses the free energy in terms of a summation over Matsubara 
frequencies as follows: 

F(p e , pi) = -T In + r dep(e)T £ In ^ 

1 — Oi J —oo 



p 

+ pp e -(T\n-^ + T\n(l + e^)+T ]T In 

V 1 + Pi n=-oo 



Pi){iuj n + p - E n - e 
1 ~ Pi 



(iu n + p)G n 



Pi- (31) 



Similarly, we can evaluate the free energy in the same fashion as Falicov and Kimbali did 



/oo 
dep mt {e) In 
-oo 



1 



-/?(e-/x) 



+ T ^ In pi + (1 - Pl ) ln(l - p^] + pp e , (32) 



where p mt (e) = y 4(1 — p%) — e 2 / (2tc) is the interacting density of states for the Bethe lattice. 
We find that both forms ( |3"TD and ( |3"2"D are numerically equal to each other, but are unable to 
show this result analytically. Since the interacting density of states is known for the Bethe 
lattice, we use the computationally simpler form in Eq. ( p2[ ) in our calculations. For the 
hypercubic lattice evaluated in Section III, we employ Eq. (|3~1~D in the free-energy analysis. 

The numerical minimization proceeds in four phases: (i) First a coarse grid is established 
for p\ and p'[ and the free energy is minimized over this grid [the electron fillings are deter- 
mined by the constraints that the chemical potential is the same in region 1 and region 2 
and that p e = ap' e + (1 — a)p", with a already determined from p, = ap\ + (1 — az)p"]] (ii) 
The filling p" is fixed at it's coarse-grid minimal value, and p\ is varied on a finer grid to 
determine the new minimum; (iii) p\ is fixed at the new minimum and p" is now varied on 
a fine grid to yield a new minimal p"; (iv) p\ and p" are varied together on the same fine 
grid to determine the final minimization of the Maxwell construction. We found that the 
minimal values of p\ and p" rarely changed in step (iv) confirming the convergence of this 
method. 



10 



We plot our results in Fig. 0. The first case considered in Fig. 0(a) is the case of 
relative half filling p e = (1 — Pi)/2. In this case the chemical potential is always at zero, 
and the relative electron filling remains unchanged for all The solid line is the first-order 
transition line and the dotted line is the spinodal-decomposition temperature. The horizontal 
distance between the solid lines at a fixed temperature is a measure of the order parameter 
p\ — p". Notice how the first-order transition temperature is always close to the spinodal- 
decomposition temperature, but that the difference becomes largest at concentrations close 
to zero and one. Note further how the two curves meet at the maximum (as they must) where 
the first-order transition disappears and becomes a second-order transition at a classical 
critical point. In Fig. 0(b) we plot the phase diagram for the case with p^ = 0.65 and 
p e = (1 — p^ I '4. In this case the chemical potential changes as a function of temperature, 
and as T — > and the system is phase separating into regions with ion densities close to 
zero and one, we find that the chemical potential will lie outside of the bandwidth of the 
interacting density of states as p'[ — > 1 because the bandwidth (4a/1 — p^) becomes narrowed 
to zero. In that case, the electron density approaches zero exponentially fast, which is why 
the spinodal-decomposition temperature approaches zero so rapidly in that regime. For this 
reason, we find that the relative electron filling is not a constant in this phase diagram, as 
it approaches zero exponentially fast near p" = 1 and it is somewhat larger than quarter 
filled near p\ = 0. The two phase diagrams in Fig. 0(a) and (b) look similar in the low- 
density regime, however. This may imply that there is an analogous scaling regime for the 
first-order T c , but it does not look like there would be a universal curve for the region close 
to p'[ = 1. The numerical effort required to perform the free-energy analysis is significant, 
so a thorough analysis of any possible scaling forms for T c was not performed. 

III. RESULTS FOR THE HYPERCUBIC LATTICE 

The formalism for the hypercubic lattice is essentially unchanged from the Bethe lattice. 
The main differences are that the integrals cannot be performed analytically anymore, which 
requires results to be worked out numerically, and requires more computational effort. The 
basic framework in Eqs. (0) — (|]) is identical as before, except now the noninteracting density 
of states is a Gaussian for the hypercubic lattice. The integral for the local Green's function 
is no longer elementary, and so one needs to solve the problem iteratively as was done 
previously for the Falicov-Kimball model: (i) first the self energy is set equal to zero; (ii) 
next the local Green's function is determined from Eq. ([5]); (iii) then the self energy is 
determined from Eq. (|8|); (iv) then steps (ii) and (iii) are repeated until the equations 
converge. We can compare the results of the interacting Green's function to the form found 
before for the Bethe lattice, by approximating the interacting density of states in the same 
fashion as before: we narrow the Gaussian by the factor a/1 — and have a total weight of 
1 — pi in the density of states. When we compare the Green's function along the imaginary 
axis at half filling in Fig. |||(a) we find that that approximation works well at high energies, 
but begins to fail near zero frequency (we chose p e = 1/6, pi = 2/3, and T = 0.1). The 
solid line is the exact result and the dotted line is the approximate (band-narrowed) result. 
The infinite-f/ Green's function on the hypercubic lattice is more complicated than on the 
Bethe lattice and the simple form that describes it for the Bethe lattice no longer holds. 
This is the main reason why the hypercubic lattice is more complicated to deal with than 
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the Bethe lattice. To see this more fully, we examine the interacting density of states in 
Fig. ^|(b). The interacting density of states is determined by solving the same equations for 
the Green's function, but this time on the real axis, rather than the imaginary axis. Notice 
how the band-narrowed form y/1 — pi exp(— e 2 /[l — piD/y/n (dotted line) is a reasonable 
approximation to the interacting density of states (solid line) but that it is too narrow, and 
it overestimates the peak height. 

Since we do not know an analytic form for the interacting density of states, we cannot 
perform the same kind of analysis that we did before at zero temperature to see if the system 
is phase separated. But we can examine the finite-temperature phase diagrams in the same 
manner. The susceptibility diverges whenever Eq. fl2~T|) is satisfied. The bare susceptibility 
now takes the form 

1 r°° I r°° I 

x° n (X) = -- P (y)— / dzp(z) 



2n J-oo ito n + p — S n — y J-oo iu n + p — £„ — yX — z\f\ — X 2 

(33) 

where X = lim^oo Yfj=i cos(kj) for the ordering wavevector k. The bare susceptibility 
continues to assume a simple form for the same three cases: (i) X = the local susceptibility 
where x°(0) = — G 2 n \ (ii) X = — 1 the "antiferromagnetic" susceptibility, where Xn(~ 1) = 
—G n /{iuj n + p — £ n ); and (iii) X = 1 the uniform susceptibility, where 

X° n (l) = = 2[1 - (iu n + p- £ n )Gy. (34) 

If we try to approximate the transition temperature by substituting in the approximate 
form we have for G n derived by assuming the interacting density of states has the same 
shape, but is band narrowed, we find that the T c 's generated are not accurate at all. Hence, 
the simple band-narrowing approximation works reasonably well for the Green's function, 
but is poor for the susceptibility. 

Instead, we simply solve for the transition temperatures numerically. We find in every 
case that we examined that the transition temperature is always highest for X = 1 and 
vanishes for all X < 0. This is shown in Fig. |4] for the cases of p e = 1/6, 1/12, 1/24, 
1/48, 1/96, 1/192, 1/384, and 1/768 and pi = 2/3, which ranges from relative half filling to 
the low-density regime. We plot T C (X) and see that the system always favors the uniform 
charge- density wave, signifying that the system wants to phase separate. We calculate the 
spinodal-decomposition temperature for phase separation by finding the temperature at 
which the uniform susceptibility diverges as a function of p e and pi. These temperatures are 
plotted in Fig. |5](a). This plot looks similar to what we found for the Bethe lattice before, so 
we try the same scaling form in Fig. |^(b), plotting TJ (pj-y/1 — Pi) versus p e / (1 — p^. Once 
again we see a data collapse, but the spread in the T c 's is somewhat larger than that seen 
in the Bethe lattice. 

Finally, we calculate the full phase diagram for the case of relative half filling p e = 
(1 — pi)/ 2 where p — in Fig. |6|. The form of this result is similar to what was seen in 
the Bethe lattice. The first-order transition temperature and the spinodal-decomposition 
temperature meet at the peak of the curve where the first-order transition becomes second 
order. We did not perform a free-energy calculation at other relative fillings here, because 
the numerical solution was significantly more difficult due to the fact that we needed to use 
Eq. ( |3~T1) rather than the computationally simpler Eq. (|32D . 
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IV. CONCLUSIONS AND DISCUSSION 



We have provided an exact solution to the spinless Falicov-Kimball model in the strongly 
correlated limit of U = oo. We only considered the less than unit density cases < p e < 1 — 
Pi, because they all satisfy a similar functional form. On the Bethe lattice we found that the 
system always phase separated at T = to states where the electrons all moved to one part of 
the lattice, and the ions moved to the other part. The spinodal-decomposition temperature 
for segregation solved a simple transcendental equation, which we showed collapsed onto a 
scaling curve. In addition, we solved for the first-order transition temperature for a select 
number of cases and discovered that the first-order transition usually occurred quite close 
to the spinodal-decomposition temperature. On the hypercubic lattice, we found similar 
results, but had to carry the analysis out numerically for all cases considered. We were 
able to explicitly show that phase separation precluded incommensurate (or commensurate) 
charge- density- wave order for the hypercubic lattice. 

These results show that when the screened Coulomb interaction is large, or in the alloy 
picture, when the A ions are extremely different from the B ions, then the system will segre- 
gate at low temperatures. This proves the segregation principle for the infinite-dimensional 
limit, and leads us to believe that it holds for all dimensions (since it has also been demon- 
strated in one dimension). Future problems to be investigated include a study of the case 
Pe + Pi — 1) as well as finite Coulomb interaction U. In the unit-density case, we expect 
charge- density- wave order to be more prevalent, perhaps precluding the segregated phase for 
all U. When the strength of the Coulomb interaction is reduced, we expect the segregated 
phase to gradually disappear and be taken over by other phase-separated or charge-density- 
wave ordered phases. 

It is possible that the phenomena described here incorporates the relevant physics of 
the charge-stripe phases in the cuprate materials: that the stripes occurred because of the 
strong propensity towards phase separation in the strongly correlated limit. The analogy 
would stem from considering the down spin electrons of the Hubbard model to be frozen 
in a particular configuration, and then examine how the mobile up spin electrons react 
to the down spins. The quantum fluctuations of the Hubbard model are replaced by the 
thermal fluctuations of the Falicov-Kimball model, and it can be viewed as a simplifying 
approximation to the charge dynamics of the strongly correlated Hubbard model, but not 
incorporating the spin dynamics. In this postulated by Emery and Kivelson,H the 

stripes would form from a balance between the desire for the system to phase separate, 
and the long-range Coulomb interaction, which would prevent the electrons from completely 
separating from the ions. There is evidence for alternative points of view, however. White 
and collaborators^ have shown that the Hubbard model on a ladder displays charge-stripe 
order even without the long-range Coulomb interaction. This order arises from the correla- 
tion of the spins and the holes and a desire to reduce the frustration induced by the hole 
motion. In their picture, the stripe ordering arises completely from a model that includes no 
long-range forces. Nevertheless, it is our belief that the phase separation exhibited here will 
be an important element of a complete description of the charge-stripe order in the cuprates 
and nickelates, because it must occur if U becomes large enough. 
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APPENDIX A: DERIVATION OF THE TRANSCENDENTAL EQUATIONS FOR 

T c ON THE BETHE LATTICE 



The derivation of Eqs. (|26| ) and (p9|) involve summations over Matsubara frequencies, 
which are performed with the help of the identity 



/3x 2 ~ 1 
tanh — = — > , 



(Al) 



for any real number x. Using Eq. Q25D and the change of variables e — > e/ y/1 — p^ yields 



Gr 



1 - l>, I —Pb{ t== 

oo e y/ 1 - pi 



iu n + p 

Employing Eq. ( |Al|) then shows that 



1 



1 



iu n + p — e iu n + p 



E 



2^1 - Pi f°° de . e 

—Pb{ n = 

oo e VI - pi 



„^oo + P P 

Eq. ( p6|) then follows from the trigonometric identity 



tanh 



(5{p 



tanh 



(3p 



tanh 



^~ e) tanlA- 

2 2 cosh 2 &±(1 -tanh 4^ tanh ' 



tanh^ 



(A2) 



(A3) 



(A4) 



The derivation of Eq. ( P§| ) is more involved, but proceeds along the same lines. Using the 
integral representation in Eq. (|25|) for G ra , performing a decomposition into simple fractions, 
and then using the identity in Eq. ( |A1| ) produces both 



E 



(zw n + /i)G n 







{iuj n + p) 2 - 4(1 - Pi ) 4 





— J\ — Pi tanh 



/?(// - 2VT^") 



-yl - A tanh 



/3(p + 2vT^ 



-oo 2-e 

oo 



cfe 



Pb{z) _ PBje) 
2 + e 2-e 



oo 2 + e 

tanh ^-6yT^)_ ; (A5) 



and 

oo 

E 



2(1 — Pi 







(iu n + p) 2 - 4(1 - P i) 



1 - Pi 



tanh 



^(/x-2vT^- 



tanh 



{3{p + 2^T=p- 



(A6) 
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Now we use the fact that the integrals for the noninteracting Green's function are trivial 

de^ = r de^ = 1, (A7) 
) 2-e J-oo 2 + e v ' 



and subtract Eq. ( JAIS] ) from Eq. ( |A6| ) to obtain 

t ^n + ,)G -2(1- p ) = Pr Y - 7 f* e tanh /WWl-/f (Ag) 

Eq. ( p9|) then follows from the change of variables e = 2COS6 1 . 
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FIGURES 



FIG. 1. Second-order transition temperature on the Bethe lattice (corresponding to spinodal 
decomposition), (a) Transition temperature plotted as a function of electron filling, (b) Transition 
temperature plotted on a scaling curve as a function of relative electron filling. 

FIG. 2. Transition temperature to phase separation on the Bethe lattice, (a) The case of 
relative half- filling (p e = [1 — Pi]/2). The solid line is the first-order transition temperature and 
the dotted line is the spinodal decomposition temperature. Notice how these two curves meet at 
the maximum where the first-order transition becomes second order, (b) The case near relative 
quarter filling (as described in the text). Notice how the shape of the curve differs from (a) near 
Pi = 1. This is because the electron filling becomes exponentially small once the chemical potential 
lies outside of the interacting bandwidth. 

FIG. 3. Comparison of the band-narrowed approximation with the exact result on the hyper- 
cubic lattice, (a) The Green's function at relative half filling on the imaginary axis. The solid 
line is the exact result and the dotted line is the approximation. The parameters are p e = 1/6, 
Pi = 2/3, and T = 0.1. (b) The interacting density of states (which is temperature-independent) 
for the exact (solid line) and approximate (dotted line) cases with p e = 1/6 and pi = 2/3 

FIG. 4. Plot of the transition temperature versus the ordering wave vector X(k) for the case 
pi = 2/3 and values of p e ranging from relative half filling (p e = 1/6) top curve, to the low-density 
regime (p e = 1/784) bottom curve, with the density reduced by a factor of 2 for each case. 

FIG. 5. Second-order transition temperature on the hypercubic lattice (corresponding to spin- 
odal decomposition), (a) Transition temperature plotted as a function of electron filling, (b) 
Transition temperature plotted on a scaling curve as a function of relative electron filling. 

FIG. 6. Transition temperature to phase separation on the hypercubic lattice for the case of 
relative half-filling (p e = [1 — Pi]/2). The solid line is the first-order transition temperature and 
the dotted line is the spinodal decomposition temperature. Notice how these two curves meet at 
the maximum where the first-order transition becomes second order. 
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